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^jl ' In the framework of the variational principle the canonical variables describing magnetohydrody- 

f^ , namic (MHD) flows of general type (i.e., with spatially varying entropy and nonzero values of all 

f^ ' topological invariants) are introduced. It is shown that the velocity representation of the Clebsch 

rvj ' type following from the variational principle with constraints is equivalent to that resulting from 
the generalization of the Weber transformation performed in the paper for the case of arbitrary 
MHD flows. Using such complete velocity representation enables us not only to describe the gen- 
eral type flows in terms of single-valued functions, but also to solve the intriguing problem of the 

' ] "missing" MHD integrals of motion. The set of hitherto known MHD local invariants and integrals 

\l . of motion appears to be incomplete: for the vanishing magnetic field it does not reduce to the set 

^~~^ • of the conventional hydrodynamic invariants. And if the analogs of the vorticity and helicity were 

^.^ I discussed earlier for the particular cases, the analog of Ertel invariant has been so far unknown. It 

(— I , is shown that all "missing" invariants are expressed in terms of the decomposition of the velocity 

^.,^' representation into the "hydrodynamic" and "magnetic" parts. In spite of the nonunique character 

i^j ' of such representation it is shown that there exists a natural restriction of the gauge transformations 

I , set allowing one to make the invariants gauge independent. It is found that on the basis of the new 

^ ■ invariants introduced a wide set of high-order invariants can be constructed. The new invariants 

(^-H ' are relevant both for the deeper insight into the problem of the topological structure of the MHD 

^ , flows as a whole and for the examination of the stability problems. The additional advantage of 

^ ■ the proposed approach is that it enables one to deal with discontinuous flows, including all types of 

' ^ ' possible breaks. 

(-H , PACS numbers: 04.20Fy, 47.10.+g, 47.65 

£?■ 

! I. INTRODUCTION. 

>■ 

^: 

lO I It is well-known that description of the solid media flows in terms of the canonical (Hamiltonian) variables is 

^^ i very useful and effective, see for instance [l|, Q. On the basis of the Hamiltonian variables it is possible to deal 

^^ ' with all nonlinear processes in unified terms which are independent of the specific problem related to the media 

~l , under investigation. For instance, all variants of the perturbation theory are expressed in terms of different order 

^-v ■ nonlinear vertices, which along with the linear dispersion relation contain the specific information relating to the 

"* — , concrete system, cf. Refs. [^|j{. In studying nonlinear stability problems the conventional Hamiltonian approach 

fx ' based upon the corresponding variational principle allows one to use the Hamiltonian along with other integrals of 

•^H ] motion (momentum, number of quasi-particles, topological invariants) in order to construct the relevant Lyapunov 

^j^' functional, cf. Refs. lMlyi0llily|- Therefore, it is important to address the problem of introducing the canonical 

r^ I variables and corresponding variational principle for the general type MHD flows (i. e., rotational, non-barotropic 

rS |. and including all types of breaks that are possible for MHD) and obtaining the complete set of the local invariants, 

L| ' see definition and discussions in original papers |lOl llllll3.ll3l | and in the recent review [l|. As for the first item, the 

• i-H , example of the variational principle describing all possible breaks is given in the recent work [lj| . 

^^ ■ Here in the framework of some modification of the variational principle of the latter work we examine the problem 

Jh ' of the MHD invariants. Note that the set of invariants for MHD discussed in the literature has been so far incomplete. 

. - - 1 The statement becomes apparent when it is considered that for the vanishing magnetic field this set has to be reduced 

to that of the conventional hydrodynamic invariants. But this limit transition does not reproduce Ertel, vorticity and 

helicity invariants existing for the hydrodynamic flows. 

Despite the fact that for the dissipation- free MHD flows there exist additional topological invariants, namely, 
magnetic helicity and cross-helicity, introduced in the papers, [ia,ll3ll3j the analogs of the vorticity and helicity 
invariants have not been discussed with necessary completeness thus far, cf., for instance, the recent review []|. The 
related quantities were mentioned for the specific cases of symmetric flows in the works [13, 1191 l20l | , the vorticity and 
helicity invariants for the incompressible flows have been obtained recently in Refs. |3, Q- I3ut analog of the Ertel 
invariant have not been presented so far (see the short communication in J2l|). The problem of obtaining the analogs 
of the hydrodynamic invariants consists in the non-potential character of the Lorentz force. Therefore, the vorticity 



and hclicity of the total velocity field v are not conserved along with the Ertel invariant construction, p~^u} ■ Vs. 
Nevertheless, corresponding generalizations have to exist, which becomes evident from the simple consideration. 
Namely, let us consider the well known set of invariants for dissipation-free MHD flows (energy, magnetic and cross 
helicity). Setting the magnetic field zero we arrive at zero values of the magnetic and cross helicity invariants, but 
do not get Ertel invariant (and hydrodynamic vorticity and helicity for the barotropic flows). This fact indicates 
incompleteness of the MHD invariants set. Evidently, there have to exist MHD analog of the Ertel invariant passing 
on to the hydrodynamic Ertel invariant for the vanishing magnetic field. Below we derive the MHD generalization 
for the Ertel invariant and show that the generalized vorticity and hclicity invariants also exist for the compressible 
barotropic MHD flows. The possibility of obtaining these invariants is based upon the velocity decomposition in 
the two parts, "hydrodynamic" and "magnetic" . The latter vanishes with the magnetic field vanishing and can be 
presented in the form of the vector product of the magnetic field and the canonically conjugate momentum and was 
first introduced in the paper [221 . In spite of the artificial character of the velocity field decomposition at first sight, we 
show that the decomposition naturally follows both from the least action principle in the canonical variables and from 
the partial integration of the Euler equations of motion (generalized Weber transformation, cf. Ref. |23,l23)- ^^^ ^^^ 
incompressible flows the latter was presented in the papers [a, 1^ . Note that the "hydrodynamic" part of the vel ocity 
is of the Clebsch type but involves vector potentials instead of the scalar ones, see discussion in Refs. |ll IliL I25ll26l| . 
The use of the vector Clebsch potentials allows one to deal with the flows possessing nontrivial topology, contrary to 
the restriction to the scalar potentials. If the latter are single-valued then the helicity vanishes identically. 

The outline of the paper is as follows. In section^ we briefly discuss the appropriate variational principle, intro- 
ducing the Clebsch type velocity representation by means of constraints and defining the canonical variables. The 
velocity representation thus obtained yields the necessary decomposition. In section ITTTl we develop generalization of 
the Weber transformation and show that it leads to the velocity representation, which is equivalent to that following 
from the variational principle under discussion. In section Hvl we examine the MHD integrals of motion, introducing 
"missing" MHD invariants, and discuss their transformation properties relating to the gauge change. We show that 
there exist natural gauges under which the additional basic invariants become unambiguous, specifically that with a 
vanishing initial value of the magnetic part of the velocity representation. In section we make some conclusions 
and formulate problems to be solved later. 

II. VARIATIONAL PRINCIPLE AND CANONICAL VARIABLES. 



Let us briefly describe the variational principle and subsidiary variables describing dissipation- free MHD. Starting 
with the standard Lagrangian density 

^ = Py -p£(p>s) + -^> (1) 

where p, s and e(p, s) are the fluid density, entropy and internal energy, respectively, H denotes the magnetic field, 
we have to include the constraint terms in the action A. Then the action can be presented as 

A^ f dtL' , ^' = / drC , jC' ^jC + Cc, (2) 

where C^ is the part of the Lagrangian density respective for the constraints, 

/dA ,. _.\ H-curlA 



Cc = pDifi + XDfi + cfDs - M • — V X curlA + VA . (3) 

\ at J An 

Here Z) = 9t + v • V is the substantial (material) derivative and A is the vector potential. [32| Including the terms with 
A and H = curlA into Cc allows us to obtain the dynamic equation for the vector potential in the gauge invariant 
form (see Eq. (|l()|l below) and to introduce relation H = curlA strictly into the variational principle. 

Supposing first that all variables introduced (including velocity) are independent, we obtain the following set of 
variational equations 

Sip =^ dtp + div(pv) = 0, (4) 

6p=^ Dip^w-v^/2, (5) 



SX=^ D/x = , (6) 

Sf^ni =^ dtXm + div(A„v) = 0, (7) 

Sa =^ Ds = 0, (8) 

,5s => dtcr + div(crv) = -pT, (9) 

SM=^ dtA = v X curlA - VA, (10) 

^A =^ a^M = ^^^^^ + curlfv X Ml. (11) 

An 

(5H ^ H = curlA, (12) 

SA =^ divM = 0, (13) 

where w and T are the enthalpy density and temperature. (3S| 

Note that in this section we assume the velocity field to be independent of other variables. Therefore, variation 
with respect to v results in the velocity representation 

Sv =^ /9v = -p\7ip - Xm\/fi,n - crVs - H X M. (14) 

It is convenient to rewrite it in a shortened form that emphasizes its structure. Bearing in mind that the velocity 
potential f, the vector Lagrange markers fi, entropy s and the vector potential A can be treated as generalized 
coordinates, one can see that p, A, a and subsidiary field — M are conjugated momenta, respectively. Let[34| 

Q^{Q,A), Q^{ip,fi,s), V = SA/S{dtQ), P=(p,A,a), r=iP,-M). (15) 

Then the velocity representation takes the transparent form 

P H X M 

v = vo(7',VQ), vo=v/,+v7(/, v/i = VQ, va/ = . (16) 

P P 

Here subindcxes h and M correspond to the "hydrodynamic" and "magnetic" parts of the velocity field. The suffix 
zero on v underlines the fact that Vq is supposed to be the dependent variable, it is expressed in terms of the canonical 
variables by means of the representation found. The hydrodynamic part, v^, corresponds to the generalized Clebsch 
representation, cf. |l4l l25l l26l|. and the magnetic part, va/, coincides with the traditional term if we replace the 
divergence-free field M by curlS. This term was first introduced by Zakharov and Kuznetsov, cf. Ref. |22|. But they 
used the incomplete form for the hydrodynamic part of the velocity, restricting it to the scalar Clebsch variables. This 
reduced form[33 evidently corresponds to the flows with zero-valued generalized helicity (or, for H ^ 0, it results 
in the hydrodynamic helicity vanishing) if the scalar Clebsch variables are single-valued. Besides, the above velocity 
representation involves the entropy term, —a\7S/ p. The latter is essential for the discontinuous flows with any types 
of the dynamically allowable breaks, including shocks, cf. Refs. [ij, lU [2J, |2g . Note that Lagrange markers /Lt are 
continuous crossing the break surface, contrary to the entropy. Therefore, the entropy term can be omitted for the 
continuous flows when the entropy can be considered as a continuous function depending on the Lagrange markers. 

From the velocity representation Eq. (|16|l and the equations of motion ^ - Hll|l it strictly follows that the velocity 
field V = Vq satisfies Euler equation with the magnetic force taken into account. Namely, providing differentiation we 
have 

„ „ curlH X H , , 

pDvo = -Vp + — , (17) 

where p is the fluid pressure. 



Canonical variables. 



The variational principle can be easily reformulated in the Hamiltonian form. Excluding the magnetic and velocity 
fields by means of Eqs. (|12|) . H16|) we arrive at the following Hamiltonian density 

n = H{V, VQ) ^VdtQ-C = p^ + peip, s) + ^ ^ ' - M • VA. (18) 

Z OTT 

Equations of motion (0J - (|11|) can now be expressed in the canonical form 

dtQ^Sn/Sr, dtV^-5H/5Q, Q := ((p, /x, s; A) , P = (p, A,a; -M) ; (19) 

Eq. (|12(l serves as a definition of the magnetic field, and the divergence-free condition for the subsidiary field M, 
Eq. H13(l , follows from the variation of the action 

A^ I dt f dr{VdtQ-n) (20) 

with respect to A. Note that it is possible to put A = 0. Under this assumption the divergence-free condition for the 
field M vanishes. But from Eq. (|ll|l it follows that divM is a conserved quantity, 9fdivM = 0. Therefore, supposing 
that divM — holds for some initial moment we arrive at the conclusion that this is valid for the arbitrary moment. 
Nevertheless, it proves convenient to deal with A 7^ that makes it possible to use different gauge conditions for the 
vector potential. 

The above variational principle results in the set of dynamic equations. From the latter follow the conventional 
MHD equations, Q, ((SJ, H17|l and the equation for the magnetic field, which follows from Eq. (|10() after taking curl, 

dtU = curl[v X H]. (21) 

On the contrary, if at some initial moment, f = i, we have the conventional MHD fields p, s, v and H, then we 
can find the initial subsidiary fields ip, fi, A, d, A, M and A, satisfying Eqs. H12(l - H14() . This can be done to within 
the gauge transformations (the latter do not change both the velocity and the magnetic field) due to the fact that 
the subsidiary fields play a role of generalized potentials. Then, if the uniqueness conditions are satisfied both for 
the conventional MHD equations and for the set of variational equations, we are led to conclude that corresponding 
solutions coincide for all moments. In this sense we can state that these sets of equations are equivalent, cf. Ref. |j|. 

The complete representation of the velocity field in the form of the generalized Clebsch representation H16() allows, 
first, to deal with the MHD flows of general type, including all types of breaks, cf. Ref. [Ij|; second, for the zero 
magnetic field it results in the correct limit transition to the conventional hydrodynamics, cf. Refs. |25l. l2fil|: third, it 
allows obtaining the integrals and invariants of motion for the MHD fiows additional to the known ones: for instance, 
the generalized Ertcl invariant, generalized vorticity and generalized helicity, see below. The two last integrals were 
deduced for the particular case of incompressible flows in the papers J3, |^ , cf. also papers [iM llil . l2(l | where the 
vorticity and helicity analogs were obtained for the MHD flows with the specific spatial symmetry. 

Moreover, it is possible to show that representation H16f) is equivalent to that following from the Weber transfor- 
mation, cf. Refs. |23ll24l | and the recent review [l|. 

III. GENERALIZED WEBER TRANSFORMATION. 



Suppose that the fluid particles are labelled by Lagrange markers a = (ai, 02, 03). The label of the particle passing 
through point r = (xi, a;2, 2:3) at time t is then 

a = a(r,t), Da = ^ + (v • V)a == 0. (22) 

at 

The particle paths and velocities are given by the inverse function 

r==r(a,t), v = i?r(a, i) = ar/5iU,_, . (23) 

Let the initial position of the particle labelled a is X, i.e., 

r(a,0)=X(a). (24) 



A natural choice of the labels would be X(a) = a; however it is convenient to retain the extra freedom represented 
by the "rearrangement function" X(a). 

We seek to transform the equation of motion l|17() to an integrable form, by generalizing the argument of Weber 
|23| (see, for example, Refs. [23, [l|, and Q. It is convenient to represent Eq. l(T7|l as 

Dv = -Ww + TVs+jxh, (25) 

where h = H/p and the vector j is defined according to 

j ^ =^. m 

being proportional to the current density. Multiplying Eq. (I25|l by dxk/dai we have 

dxk dw dxk ds dxu , .. , , dxk ,„„. 

aui oxk oai oxk oat oai 

The l.h.s. can be represented as 

iDvk)'-P^^D{J-p-)-^i.V2l (28) 

oai \ Otti ) oai 

where we have taken into account that operator D = 9/(9i|a=con5t and therefore Dxk — Vk and D commutes with 
derivative d/dui. Eq. (|27|l now takes the form 

i^(;..^)=^(.V2-H+T|^ + [Jxh].^. (29) 

\ oa.i J Oa.i oai oai 

It is convenient to transform the last term by means of the dynamic equation for the subsidiary field m = M/p 
(compare Eq. lfTT|l ) 

i:»m= (m- V)v+j/p. (30) 

Then we can transform the last term in the r.h.s. of Eq. H29() to the form of the substantial derivative, see Appendix, 

[jxh],^=i^f[mxH],^). (31) 

oa^ \ da, ) 

Analogously, the first two terms in the r.h.s. of Eq. 1)29(1 can be presented as substantial derivatives by means of 
introducing subsidiary functions (p and a, which satisfy equations (compare Eqs. ©, lO) 

D (-\ = -T, (32) 



Dip^w-v^/2. (33) 

Then 



dat dai \p J \dai p J da, \ dai J 

where we have taken into account that Ds = along with D{ds/dai) = 0. Therefore, we can present the Euler 
equation 129|) in the integrable form 



Integration leads to the relation 



dxk dip c)s a dxk . 

^fcTj— = -75 15 [Hxnifc-— +6„ (36) 

dai dai dai p dai 



Here b = b(a) does not depend on time explicitly, Dh = 0, presenting the vector constant of integration. Multiplying 
this relation by dai/dxj allows reverting from Lagrangian (a, i), to the Eulerian, (r, f), variables, 

V = -V</7 + bfcVflfe - -Vs - h X M. (37) 

P 

This representation obviously coincides with the above discussed Clebsch representation if one identifies b with 
— A/p and a with fi. Moreover, this proves equivalence of description of the general-type magnetohydrodynamic 
flows in terms of canonical variables introduced and the conventional description in Lagrange or Eulcr variables. 
The equations of motion for the generalized coordinates and momenta follow now from definitions of the subsidiary 
variables a, m = M/p, a, f and b. 

Emphasize that the vector field M = pm introduced by Eq. H30() satisfies the integral relation 



ft / M • dS = / j • dS, (38) 

where S is some oriented area moving with the fluid. This fact was first indicated in Ref. |^ for the incompressible 
flows. Now we see that it holds true for the general case. The proof of this statement is given in Appendix. Expressing 
M = curls and making use of the Stokes theorem we conclude that time derivative of the vector S circulation over 
the closed frozen-in contour 9S is proportional to the current (recall, j = (47r)~^curlH and differs from the current 
density by the constant multiplier) intersecting the surface defined by this contour, 

dt j Sdl^ f jdS^ (47r)-i [ Hdl (39) 

that highlights the physical meaning of the subsidiary field S usually introduced for the canonical description of 
MHD flows. Underline that this identity strictly follows from the dynamic equation for the subsidiary field S and is 
insensitive to the compressibility. 

The vector constant of integration, b, may be expressed in terms of the initial conditions, 

bt = v/fc(a)— h -^ h co^— , 

oai OQi oa.i 



l^fc(a) = T4.(a) + [hoxMo]fc, ^^(a) = Wfc(a,0) , (^^^ 



Lpa = (^(a, 0), Co = ( - 

\P/ t=o 

ho = ho(a) = h(x(a, 0), 0) - h(X(a), 0), Mo = Mo(a) = M(x(a, 0), 0) = M(X(a), 0). 
Under special conditions, namely, for 

X(a)=a, r(a,0)=a, a(r, 0) = r, (41) 



from Eq. igUJ) it follows 



Adopting zero initial conditions. 



we obtain 



6,:=F.(a) + ^ + co|^. (42) 



Mo-0, (^0 = 0, ao = 0, (43) 



b = V(a)=v(a,0)=vo(a)=v(a,0), (44) 



where symbol ~ indicates that we are dealing with the velocity field in the Lagrange description, i.e., v(a, t) denotes 
the velocity of the fiuid particle with label a at time t. Evidently, v(a, f) = v(r,t), where a and r are linked by 
relations H22|l and (|23|l for the specific choice given by Eqs. (|41|l . (|43|l . Then the velocity representation takes the 
particular form 

v = v,, -[hxM], V,, = -V(/5 + uofcVafc--Vs, vo(r) == v(r,0), a(r, 0) == r. (45) 

P 

It differs from that presented in Ref. [J by involving the entropy term. Note that existence of this term allows one 
to describe the general-type MHD flows (and hydrodynamic flows under condition H = 0) with arbitrary possible 
discontinuities, including shocks, slides and rotational breaks, cf. Ref. [Ij, l25l l2a |. One can omit this term for 
continuous barotropic and isentropic flows. 



IV. INTEGRALS OF MOTION. 

The conservation laws, as is well-known, follow from the specific symmetries of the action. Existence of the 
relabelling transformations group (first discussed by Salmon in Ref. |33) of the Lagrange markers, fi, leads to the 
integrals of motion that are additional to the energy, the fluid momentum and mass conservation. These additional 
integrals are expressed in terms of the Lagrange description of the motion, i.e., in terms of the Lagrange markers, etc. 
Therefore, as a rule, they arc gauge-dependent. The frozcn-in character of the magnetic field results in the specific 
topological integrals of motion, namely, magnetic helicity and cross-helicity, first discussed in Refs. [iMIiflllSj see 
also review |l| . Corresponding densities are respectively 

hM = A- H, (46) 

and 

he = v-ii. (47) 

To clarify the following discussion relating to the additional local invariants and integrals of motion, let us briefly 
recall the known ones. As it strictly follows from the dynamic equations, the local conservation law for the magnetic 
helicity holds true for general type MHD flows 

9t/iM + divqA/ = 0, qa/ = v/i^ - H • (A • v - A) . (48) 

On the contrary, in the general case the cross-helicity is governed by equation 

dthc/dt = -div [vhc + {w- uV2)H] + rdiv(sH) 

and is not conserved. But for barotropic and isentropic flows the pressure p = p(p) and he is conserved, 

dthc + diYqc^O, qc - v/ic + (X - wV2)H, (49) 

where X = / dp/ p. 

For the general case one more conserved quantity first discovered by Gordin and Pctviashvili, cf. Ref. [3l|, is known. 
Corresponding density is 

hp = ii- Vs, (50) 

and 

dthp + divqp = 0, qp = v/ip . (51) 

The integral conservation laws are related to the local conserved quantities. For instance, integrating hp over 
arbitrary substantial volume V we obtain conserved quantity Ip, 

drhp , dtip = 0. (52) 

V 

Note that hp /p g ives us an example of the so-called local Lagrange invariants (in other words, Casimirs), cf. 
Refs. [itt llll Il3. Il3| and [llQl- By definition they obey the following equations 

9ta + V • Va == , 9*1 -t- (v • V)I == , (53) 

5tJ-h(v- V)J-(J- V)v = 0, (54) 



dtL + (v • V)L + (L • V)v + L X curlv = 0, or, equivalently, dtJj + V(v • L) - v x curlL = 0. (55) 

Here a and I denote the scalar and vector Lagrange invariants, J is the frozen-in field, and L denotes S'-type invariant 
by terminology of Ref. |l2| . related to a frozen~in surface. To these invariants it is necessary to add the density p. jSg 



Evidently, the quantity hp / p is the a- type invariant. The Lagrange markers /x and quantities \/ p supply examples 
of the vector Lagrange invariants, the magnetic field H divided by p, h = H/p, is invariant of the J- type, gradient 
of any scalar Lagrange invariant is the 5- type invariant, 

L' = Va. (56) 

There also exist other relations between different type invariants, see Refs. |l|,0j allowing one to generate new 
invariants. For instance, the scalar product of the J and L invariants results in some scalar Lagrange invariant, 
symbolically 

a' = (J-L). (57) 

The above mentioned invariant hp / p can be obtained by means of this relation if we put J = h and L = Vs. Other 
examples are represented by relations generating J- (L)- type invariants by means of two L- (J-) type invariants, 

J' = [L X L']/p, (58) 



L' = p[J X J']. (59) 

Note that integrating of the density Km over an arbitrary substantial volume does not lead to the conserved integral. 
It is easy to check up that 



satisfies 



Im = / dvhM (60) 

Jv 

dtiM = / dS (A • V - A) Hn , i?„ = H • n , (61) 

Jav 

where integration in the r.h.s. is performed over the boundary dV of the volume V , n is the outward normal and dY, 
denotes an infinitesimal area of the surface dV . It is obvious that Zj\/ will be an integral of motion if _ff„ equals zero. 
This fact allows us to conclude that 2m becomes an integral of motion if we choose the substantial volume in such a 
way that on the boundary of the initial volume, I^jt^tp, holds equality Hn\t=to — 0- The latter condition is invariant 
of the motion: if equality i7„ = is fulfilled for the initial moment, then it holds true in the future. 

Another way to make Tm invariant consists in fixing the gauge of the vector potential A so that A • v = A. Then 
the dynamic equation for A, (|1U|) . takes the form 

dtA + V(v • A) - V X curlA = 0, 

i.e., A becomes an invariant of the L- type. Under this gauge condition the quantity hj^j/p presents the scalar 
Lagrange invariant, D[hM I p) = 0. 

As for the local conservation law for the cross-helicity, Eq. H49(l , it obviously leads to the integral conserved quantity 
Tc for the barotropic fiows but with the following restriction: integration has to be performed over the specific 
substantial volume, namely such that condition -ffnUy ^ ^ (this condition is invariant of the motion) holds, 

dtic = , Ic^ I drhc , H„\gy = 0. 
Jv 

Existence of the recursive procedure allowing one to construct new invariants on the basis of the starting set of 
invariants, see Refs. P, |3l, underlines the role of the local invariants among other conserved quantities. Although 
in terms of the Lagrangian variables (such as the markers /x) there exists a wide set of invariants, see, for instance, 
Ref. [l|, the most interesting invariants are such that can be expressed in Eulerian (physical) variables and thus are 
gauge-invariant . 

Emphasize that in the conventional hydrodynamics there exists Ertel invariant as, 

OLE = /i_e//0, liE^ijJ ■ Vs, (62) 

where u) = curlv is vorticity, 

dth,E + divq_E = 0, q^ = /i_ev. Due = 0. (63) 



The corresponding integral of motion reads 



dtlE = ., 



i 



drhp, . 



(64) 



Note that DIe = holds true for an arbitrary substantial volume V. 

The Ertel invariant density has the structure of the Eq. H57() with L = Vs, J = w/p (recall that a; is a frozen-in 
field for the barotropic hydrodynamic flows). In the hydrodynamic case there also exists the helicity invariant 



/i/^ = a; • V, 
which has a topological meaning, defining knottedncss of the flow. It satisfies equation 

dthn + divqn = 0, Qff = /iffV + (x - u^/2)u;, 
and evidently results in the corresponding integral conservation law 



dflH = , for 



'^nlgy 



0, 



Ih = / drhn 
'v 



(65) 
(66) 

(67) 



For the MHD case the vector w/p is not the £rozen-in field due to the fact that magnetic force is non-potential. 
It seems evident that for the MHD case there have to exist the integrals of motion generalizing the conventional 
helicity and Ertel invariant along with vorticity integral. These invariants are to pass into the conventional ones for 
the vanishing magnetic field. The generalization for the vorticity and helicity invariants was obtained in the paper 
[3 for the particular case of the incompressible flows. In the following section it is shown that there exists MHD 
generalization for the Ertel invariant, and results of the paper [g relating to the vorticity and helicity can be extended 
to incompressible barotropic MHD flows. 



A. Generalized vorticity. 



Let us prove that the quantity w^/p, where 



u}ji = cmlvfi 



V ( - ) X VQ 



V ( — ) X \/fi„ 



V I - I X Vs 



(68) 



is the frozen-in field ("hydrodynamic" part of the vorticity) for the barotropic MHD flows. It would be a trivial 
consequence of the fact that [L x L']/p, where L, L' are Lamb type invariants, is the local invariant of the frozen-in 
type if all quantities Q and P/ p satisfy homogeneous transport equations being a- or I- type invariants (remember 
that Va and V/,„ are L-typc invariants). But if and a/ p satisfy the inhomogeneous equations of motion. Therefore, 
let us start with equation of motion for the "hydrodynamic" part of the velocity. Differentiating (|16|1 and making use 
of relations 



D{VX) = V{DX) - (Vw„0 • d„,X 



we have 



Dwh = -D\-]-VQ--- V{DQ) + -(Vw„) • d„ 
Pj P P 



TVs - V(w - t'^/2) - W/imVWm , 



or, after simple rearrangements, 

Dwh = -Vp/p + {Vra - Vhm) ' Vw^ . 

Taking the curl of this equation results in 

dfLJh = -curl(w,„9„iV,i) + [Vp X Vp]/p^ - C\Xv\[vhm'^Vm) 



(69) 



[Vp X Vp]/p^ + CViA{VmVvhm ~ Vhrn^Vm) ■ 
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The term in the square brackets is equal to VmVvhm ~ t-'hm^Vm ~ v x ijj^ and we obtain 

dtLJh = [Vp X Vp]/p2 + curl[v X w;,] . (70) 

For the barotropic flows, p = p{p), the first term in the r.h.s. becomes zero and we can see that oJhl P is the frozcn-in 
field, 

At H = 0, (jJh corresponds to the conventional hydrodynamic vorticity. 

B. Generalized Kelvin's theorem. 

The frozen~in character of the generalized vorticity allows obtaining the strict generalization of the Kelvin's theorem 
for the barotropic flows. But with some restrictions it is valid also for the non barotropic flows. Namely, circulation 
r of the hydrodynamic part of the velocity over the closed material contour C is a constant of motion if the entropy 
s is constant on this contour. 

LT = 0, T = l Wh-d\ for s|c = const. (72) 

The proof strictly follows from the velocity representation: 

T ^ f Wh ■ d\ = (f { dip ^ — —dfijn H — ds ) = <J> ( -^dp.m H — ds I . (73) 

Jc Jc\ P P J Jc\ P P J 

Differentiating T and taking into account that Dp^ = D{X„i/ p) = we obtain 

DT ^ i dsD (-)= - f Tds = for s\c = const. (74) 

Note that for barotropic flows this result strictly follows from the fact that ust/ P is the frozen-in field. Namely, for 
any J-type invariant it can be easily proved that 

Z? / dS pJ • n = 0, 

where integration is performed over the substantial surface E. Then for pJ = ujh after applying the Stokes theorem, 
we have 

D dY.u}h-n^ D f Wh-dl^O. 



C. Generalized helicity. 



Now it can be proved that generalized helicity, hn^ defined in terms of the "hydrodynamic" part of the velocity, 

Hh =u}h-Vh, (75) 

is the integral of motion for barotropic flows. Differentiating Eq. H75|) and taking into account Eqs. (|69|l . H70|) for 
barotropic flows we arrive at the local conservation law of the form (rather cumbersome calculations are given in 
Appendix) : 

dthn + divqn = 0, q^ = /i^v + (x - wV2)'^h • (76) 

In analogy with the hydrodynamic case we can conclude that the integral helicity 2h (defined by means of Eq. H67|l ') is 
the integral invariant, moving together with the fluid if the normal component of the vorticity tends to zero, LOhn = 0, 
on the surface of the corresponding substantial volume V. Note that the condition tohn = is invariant of the motion 
(due to the frozen-in character of uj^/ p) and therefore it can be related to the initial surface only. 
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D. Generalized Ertel invariant. 

Let us show that there exists strict generahzatioii of the Ertel invariant for the MHD case. For this purpose let us 
prove that without any restrictions related to the character of the flow the quantity 

hE = [uJh ■ Vs) (77) 

obeys the conservation law of the form 

dthE + divqs = 0, q^ == Hev. (78) 

Equivalcntly, the quantity aE = hs/p is transported by the fluid 

DaE = 0, aE = fiE/p, (79) 

being a- type invariant. For the barotropic flows it immediately follows from the fact that uJhl P is the frozen-in field 
if the composition rules given by Eqs. H57(l and 1)56(1 are taken into account. In order to make the proof for the non 
barotropic flows more transparent let us consider a more general situation. Let J satisfy equation of motion of the 
form 

m = (J • v)v + z, (80) 

differing from the frozen field equation (|54|) by existence of the term Z that violates homogeneity. Then, if a represents 
any scalar Lagrange invariant, we have 



^(J • Va) = DJ • Va + J • DVa = Z • Va + ( (J • V)v 1 • Va - (J • Vw„) • d^a. 

Here the two last terms cancel and we get 

£>(J ■ Va) = Z • Va if m = (J • V)v + Z and Da = 0. (81) 

For Z := these relations prove the generating rule of Eq. ((57(1 . But we can see that J • Va becomes the local Lagrange 
invariant under less restrictive condition Z • Va = 0. That is the case for the Ertel invariant: Z = [Vp x Vp]/ p^ 
is orthogonal to Vs due to the fact that the scalar product of any three thermodynamic quantities is equal to zero 
(because any thermodynamic variable in the equilibrium state is a function of two basic variables). This concludes 
the proof. 

The conserved integral quantity associated with aE is 

1e = / drhE , dtlE = 0. (82) 

Jv 

Note that by the structure Xe is not gauge-invariant in contrast to the hydrodynamic case. Let us examine its change 
under gauge transformation that results in v/j ^ vj^, vj\/ => v^^ with 

^h + ^M = V/j + VM • 



Then 



I'e -Ie = I dr Vs • (oiJi - ^h) = I dr Vs • {u^m - i^'m) 
Jv Jv 



'V JV 

But Vs • {ujaj — w'^j) = — div[Vs x (v^^ — V7\/)] and, therefore, 

I'e-Ie = - f dEn-[Vsx (v^^-vm)]. 
Jav 

Now we can proceed in the two ways. First, making use of identity Vs x X = curl(sX) — s ■ curlX we obtain 

I'e - Ie = - dS n • (curl(s(v^j - vm)) - s curl(v^j - vm)) . 

JdV 
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Here the integral of the first term vanishes (that is trivial for a closed boundary dV and assumes the necessary decrease 
of the integrand for the infinite volume V) and we have 



X'e —Xe= dT, sn ■ curl(v^j 
Jav 



vm) . (83) 



This representation immediately suggests that integral Ertel invariant becomes gauge-independent for the substantial 
volume V chosen in such a way that its boundary coincides with the entropy-constant surface, sLy = const. 
The second way is as follows. Bearing in mind that v^^ — v^/ = — [h x (M' — M)] we obtain 



I'e-Ie^ /__ dSn- [Vsx [hx (M'-M)]]. 

JdV 



(84) 



Inasmuch as both M' and M satisfy Eq. (|ll|l . their difference, M = M' — M, is governed by the homogeneous equation 

5tM = curl[v X M], 

i.e. m = M/p is frozen-in field. Then we can conclude that the vector Vs x [h x m] entering the integrand is the 
frozen-in field, as it follows from recursion relations H56() - (|58|1 . Therefore, if we adopt relation n-[Vsx [hxTn]]|g^, = 
as the initial condition, then it holds true for all moments. For instance, this relation is fulfilled if H^ — and to„ = 
at the initial moment. Evidently, these two conditions cannot be fulfilled for an arbitrary gauge. But we can restrict 
ourselves to a such subset of the initial conditions for the subsidiary field M that M|t=tQ = /H|t=tQ, where / is an 
arbitrary function. (Then divM = divMjt^tg = V/ ■ H|t=tQ is time-independent in accordance with equation Hll|l . 
and for the particular choice of / such that H|t=tQ • V/ = wc have divM ~ 0.) For these initial conditions M along 
with M' arc collincar to H at the initial moment and therefore the initial value of the scalar product n ■ [Vs x [h m]] 
is zero. Due to the frozen-in character of the quantity Vs x [h x Tn] equation n • [Vs x [h x m]] = holds true for 
the arbitrary moment. Thus we can make the conclusion that gauge dependence of the Ertel's invariant can be partly 
eliminated by appropriate choice of the initial conditions or substantial volumes. 

E. Specific gauge. 

Examination of the integrals of motion shows that they are gauge-dependent. This dependence is attributed to 
different decompositions of the velocity field into the "hydrodynamic" and "magnetic" parts, v = v/j + vm- Underline 
that there exists a wide subgroup of the gauge transformations that include transformations that change the generalized 
potentials (p, /x, A and a with no change in M. The "hydrodynamic" part of the velocity representation, v^, does 
not evidently vary under the action of transformations of this, say, "hydrodynamic" subgroup. Then the generalized 
circulation, integral helicity and Ertel integral are likewise invariant under these gauge transformations. The simplest 
way to restrict the gauge transformations by this subgroup consists in adopting zero initial conditions for the subsidiary 
field M. This choice does not restrict in any way the character of a flow, in particular, all integrals of motion can 
possess nonzero values. The more detailed discussion of the gauge dependence of the additional integrals considered 
will be presented elsewhere. 

V. CONCLUSIONS. 

The results obtained can be summarized as follows. First, the variant of introducing the canonical description of the 
MHD flows by means of the variational principle with constraints is presented. It is shown that in order to describe 
general-type MHD flows it is necessary to use in the generalized Clebsch-type representation for the fluid velocity 
field the vector Clebsch variables (the Lagrange markers and conjugate momenta) along with the entropy term (cf. 
papers J25i l26l | describing the hydrodynamic case) and the conventional magnetic term introduced first in Ref . |22l | . 
Such a complete representation allows one to deal with general-type MHD flows, including all type of breaks, see Ref. 
|14|. Second, it is proved that the generalized Weber transformation introduced leads to the velocity representation, 
which is equivalent to that introduced by means of the variational principle. Third, the existence of the generalized 
Ertel invariant for MHD flows is proven. Forth, there are generalized the vorticity and helicity invariants for the 
compressible barotropic MHD flows (first discussed for the incompressible case in cf. Q)- Fifth, the relations between 
the local and integral invariants arc discussed along with the gauge dependence of the latter. 

As a consequence of the completeness of the proposed velocity representation we get the correct limit transition 
from the MHD to the conventional hydrodynamic flows. The results obtained allow one to consider the complicated 
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MHD problems in terms of the Hamiltonian variables. The use of this approach was demonstrated for the specific 
case of incompressible flows in the series of papers [3, |^ devoted to the nonlinear stability criteria. Emphasize that 
existence of the additional invariants proved in our paper is of high importance for the stability problems. 

Note that existence of the additional basic invariants of the motion makes it actual to examine the problem of the 
complete set of independent invariants, and, respectively, the complete set of the corresponding Casimirs, cf. Ref. [l|. 
For instance, existence of the three independent basic local invariants for the non barotropic flows (s, as and h) 
immediately leads to the two denumerable sets of the monomial scalar invariants 

(™) ;'u V7\m (™) fu V7\m (1) (0) n 1 o 

ap'~[h-\/) s, a^ = (h ■ V) as, ap' ~ ap , a^^' — as , m — 0,1,2,... 

The first set was discussed in the paper [ll , and the second subset is a new one along with the "parent" Ertel invariant 
aE- Evidently, 






where / is an arbitrary function, is also the scalar Lagrange invariant. Therefore, we immediately arrive at the 
following set of integrals of motion (Casimirs) 



drpfl{a^r},{a'->} , (85) 

V ^ ' 

which is much wider than that discussed in the literature, cf. Eq. (10.23) in Ref. Q. The additional set of the scalar 
monomial Lagrange invariants can be generated by the magnetic hclicity under the specific gauge condition, A = A-v, 

a]^ = (h • V)"aM. This enables evident generalization of the integrals of motion (|85|l . 
One example of the additional J- invariants reads 

J' = [Vs X Vaij]//9. 

In turn, one can get new sets of the scalar invariants by applying operation (J' • V) to the previous scalar invariants 
and so on. Obviously, this also leads to additional Casimirs to that indicated in (|85|l . 

For the barotropic fiows the picture is analogous: the basic set of the scalar Lagrange invariants involves the 
generalized hclicity, an = huj p, and au (under the gauge condition specified above), and with J-invariants h and 

ijJh- Therefore, we obtain additional scalar invariants a)^ ~ (h ■ V)"^//, a\" , or^ ~ {p^^ujh ■ V)"a//, 5]y" ~ 
{p^^ujh ■ V)"aA/ and the conserved integrals of the form 

drpf[{a^-^},{aii\{a^p},{a^-p}) . (86) 

This set of the Casimirs generalizes that presented in cf. Ref. [j , the latter follows from (|86|l if we replace the function 

(n) 

f depending on the four sets of the monomial invariants by a function depending only on the invariants a\.j . 

Note that we can construct the following generations of the local invariants by means of the recursion relations and 
obtain Casimirs of a more sophisticated structure than that presented in H85|) . H86|) . The problem of obtaining the 
complete set of the local invariants and gauge invariance of the corresponding integral invariants is rather complicated 
and is still under examination. This questions will be discussed in detail in the forthcoming paper. 

Appendix A 

In order to prove Eq. (|31|) let us substitute j from Eq. H3U|) into expression [j x h]kdxk/dai. Then 



[j X h],^ = [Dm X H],^ - [(m • V)v x U],^^ = 
oai oui oai 

= ^D ([m x H],) - ([m x D{ph)]k + [(m • V)v x H]^) ^. 
oai oui 

Proceeding with the terms in the second brackets we obtain 

[m X D{ph)]k + [(m • V)v x H]fc = [m x h]fc • Dp + [pm x Dh]k + [(m • V)v x H]fe = 
= -[M X h]fc • divv + [M X (h • V)v]fc + [(M ■ V)v x h]fc = -[M x h],,dkVs 



(87) 
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where M = pm and the dynamic equation Dh = (h • V)v along with identity 

[M X (h ■ V)v]fc + [(M • V)v X h]fc = [M X h]kdsvs - [M x h]Avs 
are taken into account. Introducing for brevity notation 

Y = mxH = Mxh, 
we can represent the r.h.s. of Eq. H87(l as 

oai oai oai oai oai oai \ oai 

This proves Eq. l|^ . 

Let us check up the integral relation (|38|l . It is sufficient to prove the differential form, namely 

£)(M-dI]) =j-dS, (89) 

where dS is some infinitesimal oriented area moving with the fluid. It can be presented as 

dS = dli X dh, (90) 

where dli, dh are frozen-in linear elements. Thus, dlo, a = 1,2, are invariants of the J- type and satisfy equations 

D{d\a) = {d\a ■ V)V. 

Consequently, from the recursion relation Eq. H59|l it follows that pdTi is L- type invariant and hence it is governed 
by the dynamic equation 

D{pd'S) = -V(pv ■ d£) + V X curl(pdi;), 

or in the coordinates, 

D{pdl:,)^-{pdl:k)^^Vk. (91) 

Now it is easy to prove relation (|89() without any restrictions for the type of flow. Namely, 

D{M. ■ dS) = D{va ■ pdT.) = Dm ■ pdT. + m,D{pdT.,) = 
= pdH ■ (m • V)v + j • dS — mipdT,kdiVk = j • dY, . 

In order to prove the helicity conservation, Eq. H76() , let us consider some scalar quantity of the form 

where J is some frozen-in field. Then, taking into account that Eq. (|69|) for the barotropic flows can be rewritten as 

Dvh = -V(x - v'^/2) - Vhrn ■ Vwm , X = / dp/p, 

we obtain 

DY = Dwh • J + Vft • L>J = -V(x - v'^/2) ■ J. 

For J = ijJh/ P we proceed 

D{vn ■ <^hlp) = ~p-^ (V(x - t-VS) ■ ^h) = -p-Miv ((x - v''/2)un) ■ 

Then 

D{\h ■ ^h) = pD{wh ■ i^h/p) + (v/i • u}h/p)Dp = -div ((x - v^ I2)u}h) - {^h ■ t^h) divv, 

or 

9t(v/, ■ w?j) = -divq^, q,i = (x - v^/2)a;h + v (v^ • u;^) (93) 

that evidently coincides with Eq. (|76|) . 

It is noteworthy that the proof is valid for arbitrary J- type invariant if the field p3 is divergence-free and the flow 
is barotropic: 

dtipJ ■ Vh) = -divq , q = (x - v^/2)p3 + v (pJ ■ v,,) for div(pJ) = 0. (94) 

For instance, choosing J = h immediately leads to the cross-helicity invariant if one takes into account that H • v/j = 
H V. 
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